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Abstract 

We derive evolution equations satisfied by moments of parton distributions 
when the integration over the Bjorken variable is restricted to a subset (xq < 
x < 1) of the allowed kinematical range < x < 1. The corresponding 
anomalous dimensions turn out to be given by a triangular matrix which 
couples the iV-th truncated moment with all (N + K )-th truncated moments 
with integer K > 0. We show that the series of couplings to higher moments 
is convergent and can be truncated to low orders while retaining excellent 
accuracy. We give an example of application to the determination of a s from 
scaling violations. 
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The description of scaling violations of deep-inelastic structure functions is his- 
torically one of the first predictions of perturbative QCD amenable to experimen- 
tal testing, and still is one of the most accurate ways to use perturbative QCD 
for precision measurements Specifically, comparison of the theoretical predic- 
tion with the data allows a determination of the only free parameter in the QCD 
lagrangian, the strong coupling a s , as well as the extraction of the parton distri- 
butions of hadrons which, though in principle computable, are determined from 
the nonperturbative dynamics of the theory and thus must be treated as unknown 
phenomenological parameters. 

As well known, scaling violations are described by ordinary linear differential 
equations for evolution in t = In satisfied by the Mellin moments of parton 
distributions, which can be directly viewed as matrix elements of local operators. 
At the leading log level, parton distributions can be expressed directly as linear 
combinations of measurable structure functions; beyond leading log this is only 
true in specific factorization schemes, whereas in a general scheme the moments of 
structure functions are related to moments of parton distributions through Wilson 
coefficients which are calculable as perturbative expansions in a s [0. 

From a phenomenological point of view, however, dealing with moments of struc- 
ture functions is rather inconvenient, since by definition Mellin moments are ob- 
tained integrating over all values of the Bjorken variable < x < 1. Since x is 
related to the invariant energy W 2 of the virtual photon-hadron scattering process 
by W 2 = - =s , x — > is the infinite energy limit and can thus never be attained 
experimentally. All moments are thus subject to an a priori infinite uncertainty 
from this region, which can only be reduced on the basis of nonperturbative models 
and assumptions, such as the idea that the virtual photon-hadron scattering cross 
section should behave at high energy like the cross section for scattering between 
real hadrons and thus be controlled by Regge theory ||. It follows that any use of 
the evolution equations for moments requires some model-dependent input. 

A way to avoid this problem is of course well known: work in x-space and 
deal with Altarelli-Parisi equations, which give directly the evolution of parton 
distributions (rather than of their moments). Undoing the moments turns the 
ordinary differential equations in t into integro-differential equations in x and t, 
but the x integration is such that scaling violations at xq only depend on the 
values of parton distributions for all x > xq. The need for an extrapolation to the 
unmeasurable x — > region is thus at least in principle avoided B. 

In practice, however, in order to be able to solve the Altarelli-Parisi equation, 
we must parameterize the data in some specific way. This is usually done by fitting 
a well-defined functional form f(x), parametrized by some free parameters, to data 
at a fixed scale Q 2 ; or, alternatively, expanding on suitable basis of functions 0. 
In either bias is introduced in the analysis, due to the fact that the specific 

choice of functional form of the fitting function, or of the basis functions, constrains, 
for obvious reasons of smoothness, the description of data with the smallest mea- 
sured values of x. Giving a precise quantitative estimate of the possible bias thus 
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introduced can in practice be rather complicated 0. 

One may also be interested in using scaling violations for the determination of a 
moment of a parton distribution. For instance, the gluon distribution is primarily 
determined from scaling violations pj, and one may be interested directly in its 
moments because of their physical interpretation, or because the structure of evo- 
lution equations is such that some moments are better constrained by an observed 
pattern of scaling violations than the parton distribution for any individual value 
of x ||. However, when arriving at a determination of, say, the first moment of 
the polarized gluon distribution ||, it would be very useful to be able to separate 
the contribution to the moment from the measured region — which is determined 
from the data up to conventional uncertainties — from that due to an extrapolation 
which is necessarily entirely based on theoretical prejudice. 

All these problems are solved if one is able to formulate evolution equations 
directly in terms of moments restricted to the measured region. Since, as already 
mentioned, evolution at Xq requires knowledge of the parton distribution for all 
x > xq, it is clear that the scale dependence of a moment evaluated by integrating 
over the restricted range Xq < x < 1 (truncated moment, henceforth) rather than 
over the full kinematically allowed region < x < 1, is determined by the structure 
function in the same region.^ It is however not obvious that the evolution equations 
for truncated moments will take a simple form, or even that they will close upon a 
specific subset of moments: if they did not, nothing would be gained by introducing 
truncated moments over the resolution of the full Altarelli-Parisi equations. 

Here, we will derive evolution equations for truncated moments. We will show 
that the evolution of truncated moments is driven by a triangular matrix of anoma- 
lous dimensions which couples the iV-th truncated moment to all N+K-th moments, 
where if is a positive integer (but N can be any real number). The elements of 
this matrix of anomalous dimensions depend on the cutoff xq, and are calculable 
in perturbation theory as straightforward integrals of the Altarelli-Parisi splitting 
functions. Furthermore, we will prove that the series of couplings to higher mo- 
ments is convergent, so that the infinite matrix of anomalous dimensions can be 
truncated to specified accuracy. We will also see that this convergence is very fast, 
so in practice it is only necessary to deal with small (typically less than 10 x 10) 
matrices. 

In order to derive the evolution equations for truncated moments, we start from 

2 Turning the argument around, this also means that evolution equations for moments truncated 
to a generic interval xq < x < x\ do not close, because they would also require knowledge of parton 
distributions for x > x\. Hence, problems related to the extrapolation to x = 1 of data taken only 
at x < x\ cannot be solved by the methods of this paper. The extrapolation to x = 1 is however 
a much less serious problem, both because the x — > 1 limit is, unlike the x — > limit, at least 
in principle experimentally accessible, and because structure functions must vanish as x — > 1 for 
kinematic reasons, so that the extrapolation is under much better control. 
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the Altarelli-Parisi evolution equation 



at Z7T Jx y \y ) 



Here q(x, Q 2 ) is the nonsinglet quark distribution; in the singlet case one must 
consider a 2 x 2 matrix of anomalous dimensions which mixes the quark and gluon 
distributions. This introduces some trivial complications which we will not discuss 
here. The splitting function P(x) is given as a series in a s , 
henceforth we will suppress its explicit a s dependence. 

The truncated Mellin moment of the parton distribution q(x, Q 2 ) is defined as 

q N ( Xo ,Q^) = f 1 dxx N - 1 q{x 1 Q 2 ) . (2) 

JXQ 

By integrating Eq. (|I|) and inverting the order of the double integration it is easy 
to see that truncated moments satisfy the evolution equation 



4v(* ,Q 2 ) = ^f 1 f ~ dyy N - x q{y,Q*)G N (^;a s (Q 2 )) 

at Z7T Jx \ V J 



(3) 



with an evolution kernel given by a truncated moment of the splitting function 

G N (x) = f 1 dzz N - l P{z) . (4) 



If x = the kernel G^{xo/y) reduces to the usual x-independent anomalous 
dimension, Gat(O) = ^n, it can thus be taken outside the integral in Eq. (^), and 
the r.h.s. of Eq. ([3]) depends only on the iV-th moment. Henceforth we will assume 
that ReiV is large enough for 77V to be regular. If instead xo 7^ 0, because of the 
residual y dependence in the kernel Gn the evolution equation does not diagonalize. 
However, it is easy to see that the N-th truncated moment mixes only with moments 
with index M >N. 

To prove this, expand the kernel GV in Taylor series around y = 1. 

^0 nl 




where 




(6) 



Since Gn{z) is regular for all < z < 1, but in general has logarithmic singularities 
as z — > 1, due to the presence of + distributions in the splitting function P(z), 
the Taylor expansion in Eq. ([5]) has radius of convergence r = (1 — xq). However, 
since the singularities of Gjvt^o/l/) at y = xq are integrable, we can substitute the 
Taylor expansion Eq. (|||) in the r.h.s. of the evolution Eq. (||), exchange the order 
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of sum and integral, and still end up with a convergent sum. Since Gn{xq/u) is 
regular at y = 1, the expansion Eq. (|5]) contains only non-negative powers of y, so, 
after substitution of the expansion in Eq. (Q), there is no mixing between the iV-th 
truncated moment and moments with M < N, as promised. 

Because of the convergence of the Taylor expansion (||) and of the ensuing ex- 
pansion of the r.h.s. of Eq. @, we can truncate the expansion at finite order M. 
A straightforward computation then leads to 

G N (-) = lt c ( K,kxo)y K + 0\(y- , (7) 



K=0 



where 



v=K K\(p-K)\ ' 
so that the evolution equation (Q) becomes 

^-q N {xo,Q 2 ) = - E c( ^n( x o}Qn+k(x ,Q 2 ) • (9) 

at 27r ^— ' 

We are thus led to an ordinary finite system of differential equations, which can 
be solved by standard methods, provided the number of equations in the system 
equals the number of unknowns. For this to happen, we must include a decreasing 
number of terms in the series as the order of the moment increases. For example, 
if one is interested in the evolution of the iVo-th moment, and wishes to include 
M + 1 terms in the series that expresses its scale dependence, one must include in 
the series associated with the (Nq + i^)-th moment only M + 1 — K terms. This 
then gives an upper triangular matrix of coefficients. Such an approximation is 
only possible if higher moments have a decreasing influence on the evolution, so 
they may be approximated less accurately. 

It is easy to see that this is indeed the case. If the Taylor expansion Eq. (H) is 
truncated at order M, the percentage error on the r.h.s. of the evolution Eq. 
due to the truncation of the series, is equal to 



R(N, M; x , Q 2 ) = 1 £ dyy N ~ l q{ yi Q 



y . 

where the normalization is given by the exact integral 



K=0 



, (io) 



Af = jjyy N -\{y,Q 2 )G N (^ . (11) 

Now, the coupling of the N— th moment to the (N + M)— th moment is due to 
terms which are at least of order M in the Taylor expansion Eq. @ (because the 
expansion of the coefficient ck,n in Eq. (|D starts at order K). Such terms decrease 
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very rapidly because the Taylor series is convergent, and furthermore the relative 
size of the (N + M)— th moment compared to the iV-th moment decreases rapidly, 
since parton distributions fall off as a power of (1—x) as x — > 1. We will check this 
explicitly below in the particular case of the NLO evolution of the nonsinglet quark 
distribution, however we emphasize that it follows from rather general properties 
of parton distributions and their evolution. 

Once the system @ has been truncated to a finite size, it is straightforward 
to solve the evolution equations explicitly, by techniques analogous to those used 
to solve the standard coupled singlet evolution equations. At leading order, the 
evolution kernels are t-independent, so the only t dependence on the r.h.s. of the 
evolution equation is in the explicit factor of a s (Q 2 ). The solution is thus simply 
obtained by diagonalizing the matrix of coefficients, a task which is in our case is 
enormously simplified by the fact that the matrix is upper triangular. 

Let us derive the leading order solution to the evolution equation of the iV -th 
truncated moment, using Eq. (|9|), which we can rewrite in simplified notation as 

EC, (12) 

dt 2n L=N 

where Nq < K, L < Nq + M, and the matrix of coefficients is given by 

= c<™°>(*„) (L>K) , 




(L < K) . < 13 > 

A few basic properties of triangular matrices are collected in the Appendix. One of 
them, useful below, is the fact that the eigenvalues of a triangular matrix such as 
C coincide with the diagonal elements, Ckk = c\^ I K ~ K+N °\ 

Define now the rotated moments, in the basis in which C is diagonal 

N +M 

q K = R klQl , (14) 

L=No 

where 

No+M 

RklClpRpq = Ckk^kq ■ (15) 

L,P=N 

The matrix R which diagonalizes C, and its inversei? -1 , are also upper triangular, 
and can both be computed exactly by means of a simple recursion relation in terms 
of the elements of C, without having to resort to the time-consuming evaluation of 
determinants (see the Appendix). It is apparent that the rotated moments evolve 
independently, and the solution to their evolution equation is given by the familiar 
expression 



q K (xo,Q" 



a s (Q 2 )_ 



q K (x 0} Q ) , (16) 
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where 60 is the leading coefficient of the (3 function, 6 = 11/2 — n//3 for SU(3). 
The solution to the evolution equation for truncated moments is then found by 
simply rotating back, 

N +M 

q K = J2 r klQl ■ 

L=N 

At next-to-leading order, the evolution kernel for truncated moments, Gn, ac- 
quires a scale dependence though a s (Q 2 ), and can be written as 



(17) 



(o) fx \ , a s (<5 2 ) G (i) / Xq 



Gn \ — \ + 



y 



2tt 



r N 



The NLO evolution equation can then be written as 



dq 



K 



dt 



a s (Q 2 ) N ^ 



2tt 



L=N 



r (o) , 



a s (Q 2 
2tt 



' U KL 



Ql 



(19) 



where the matrices C*- -* and C*- 1 -* are given by Eq. (|13|) in terms of the coefficients 
c% N and Cfr N constructed according to Eqs. (|)-|D from the LO and NLO kernels 
respectively 

Diagonalizing the matrix of leading order anomalous dimensions with the matrix 
R of Eqs. QT3| ) and ( P~5j ) we get an evolution equation for the rotated moments q, 
of the form 



dq 



K 



dt 



a s (Q 2 ) N ^ 



2tt 



L=Nq 



r a i_ a *(Q 2 ) h 

^kkOrl H u kl 

In 



Ql 



(20) 



where 



(21) 



N +M 

Drl = RkpC P qRq L ■ 

P,Q=N 

The matrix evolution equation (PH|) can be solved with standard techniques of per- 
turbation theory The evolved (rotated) moments are expressed in terms of the 
initial condition and of the various anomalous dimensions involved as 



<k(Q 2 ) 



\<Xs(Q 2 o)~ 




Ws{Q 2 )_ 





N +M f) I 
y> L>KL A 

L=N 27T CkK — ClL + ^0 



(a s (Q 2 Q ) ~ a s (Q 2 

C LL /b 



MQ 2 ) 



a s (Q 2 



MQ 2 ) 



C'KK/bo 



Ql(Qo) 



(22) 



where b\ is the second coefficient of the QCD (3 function, b\ = 51/2 — 19n//6 for 
SU(3). The NLO solution is then found once again by just rotating back according 
to Eq. (0). 
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In order to show this formalism at work, we consider now the determination of 
a s from scaling violations of the nonsinglet structure function -F 2 NS . Our method 
allows a direct determination of the strong coupling by fitting the evolution of the 
truncated moments of the measured structure function, with a s left as the only 
free parameter, without having to introduce a parametrization of parton distribu- 
tions. A determination of a s from scaling violations of a nonsinglet moment has 
been sometimes attempted, but only in the presence of sum rules which fix the 



asymptotic normalization of the moment g, |10| . However, even the presence of this 
constraint does not obviate the problem of the uncertainty introduced by the small 
x extrapolation, which remains sizable [1], |J . 

The structure function F 2 m the DIS scheme || is simply equal to the nons- 
inglet quark distribution, 

F 2 NS (x,Q 2 ) = (F 2 p (x,Q 2 )-F 2 n (*,Q 2 



(23) 



r 

XX 2 [ii( x ^Q 2 ) +qi(x,Q 2 



i=l 



p—n 



In order to determine its evolution it is thus sufficient to transform to the DIS 
scheme the well-known NLO nonsine let splitting function P* s |TJ. We can then 
determine the evolution kernel Gn(z) Eq. (f|) analytically, and study the accuracy 
of the truncation of the expansion in Eq. (|5]) by explicitly computing the function 
R(N,M;x ,Q 2 ) defined in Eq. (|10|). For this purpose, we must use an explicit 
nonsinglet quark distribution, which we can take from any recent parton distribution 
set. We then evaluate the function R at the reference scale of the chosen parton 
set, with several typical choices of the lower limit xq of the x-range, by including 
a decreasing number of terms as the order of the moment increases, as discussed 
above. 

The results are shown in Table 1 for moments between the second and the fifth. 
It is apparent that, despite the fact that less terms are included, the accuracy 
of the determination of higher moments is actually higher: the fact that higher 
moments are largely insensitive to the lower limit of integration overwhelms the 
error introduced by the truncation of the Taylor expansion in Eq. (^). In particular, 
it is apparent that in order to reliably compute the evolution of the second moment 
it suffices to consider a four by four evolution matrix. An accurate description of 
the scaling violations of the first moment would instead require the inclusion of 
several more terms; this is a consequence of the fact that the integrand in Eq. (f|) 
decreases as z — > only if iV > 1 . The convergence to the correct result is of course 
slower when xq is larger, since the dependence on y of the kernel Eq. (§) is weaker 
when x is smaller (and indeed, there would be none in the limit Xq = 0). 

The simplest way to determine a s is to extract from the data the required set 
of moments, for instance the moments from the second to the fifth, as in Tab. 1, 
at the scale where the kinematic coverage in x is widest at large x (so all moments 
can be reliably determined), then solve the evolution equation, compare to the data 
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x 


O.Ol 


0.03 


0.1 


0.01 


0.03 


0.1 


N 


LO 


NLO 


2 


6.3 10~ 3 


3.3 KT 2 


1.5 10 1 


3.5 10~ 3 


2.7 10~ 2 


2.0 lO" 1 


3 


1.0 10~ 4 


1.7 10 3 


3.010- 2 


6.3 10~ 5 


2.8 10~ 3 


3.3 10~ 2 


4 


1.710- 6 


8.6 10~ 5 


5.1 lO" 3 


1.1 10~ 6 


6.9 lO" 5 


5.5 10~ 3 


5 


2.7 10~ 8 


4.1 10 _B 


8.3 KT 4 


1.8 HT 8 


3.3 lO" 5 


8.7 10~ 4 



Table 1: The percentage error function R(N,M;xo,Q 2 ) defined in Eq. ([TOD, com- 
puted from the LO and NLO contributions to the nonsinglet splitting function in 
the DIS scheme, with M = 5 — N, the values of iV and Xq shown, Q 2 = 2.56 GeV 2 
and nonsinglet quark distribution from the CTEQ4D parton set [pT 



for the second moment at other scales (where, for example, the coverage at large x 
might be smaller so higher moments are less accurately determined) , and perform a 
fit of a s . Our purpose here is not to perform a detailed phenomenological analisys, 
but rather to explore the viability of the method. 

We have thus simply attempted such a fit of a s by using as "data" a parametriza- 
tion of all available data on F 2 (x, Q 2 ) for proton and deuteron targets fl3fl , which 
(neglecting nuclear effects) determines the nonsinglet as _F 2 NS = 2(F 2 P — _F 2 d ); the 
moments are then simply found by numerical integration of the parametrization. 
The kinematic range is essentially limited by the availability of deuterium data: 
even with a truncation point xo — 0.1, a reliable reconstruction of the moments is 
possible only for 30 GeV 2 ~ Q 2 ~ 100 GeV 2 . Imposing that power corrections be 
negligible requires the lower cut at Q 2 > 30 GeV 2 . In fact, since power corrections 
are large at large x 0, this is important in order for the determination of the higher 
moments to be reliable. 

With all these cuts, fitting to such "data" gives a s (M z ) = 0.115. An estimate 
of the statistical error is unfortunately impossible, since a reliable determination of 



the covariance matrix for the best-fit parameters is not available |14| for the fits 



of Ref. [0.0 However, the fact that the central value is so close to the current 



global DIS average [a s (M z ) = 0.117 ± 0.002(exp.) ± 0.004(th.)] suggests that 
the statistical error is rather small. 

It thus appears worth considering an actual extraction of a s from the data using 



3 Since the parametrization of F2 is provided with an "estimated error band" one might 
hope to get a qualitative idea of the error by taking the integrals of the upper and lower curves of 
the band as estimates of the error on the moment. This procedure is however meaningless, as seen 
by noting that the error on a s could then be made arbitrarily small by increasing the number of 
values of Q 2 at which the moment is evaluated (even within a fixed range in Q 2 ). This apparently 
paradoxical result is of course due to the fact that the procedure neglects correlations between 
the values of the moment extracted from the fit at two different scales, which tend to one as the 
scales get closer. 
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this method, either by computing the moments numerically from a single set of 
data, or by first determining a global parametrization of the data with correlations 
taken into account as required. Such an extraction could presumably be improved 
by optimizing the choice of moment to be fitted. Indeed, high moments cannot 
be determined accurately from the data, due to the poor knowledge of structure 
functions at large x, while the evolution of the first moment (which is the lowest 
convergent one in the nonsinglet case) is hard to determine accurately due to the 
slower convergence of the expansion Eq. (|5|). Note that the optimal moment need 
not be integer, so one would rather expect to have an optimal range. Determinations 
of a s from different moments in this range could then be combined by properly 
taking their correlations into account. The extension of the formalism to the singlet 
sector (which is essentially straightforward) will allow both a determination of a s 
from wider data sets, and a determination of the partial moments of the gluon 
distribution, which, as we already discussed, might be of great phenomenological 
relevance. 

In summary, we have determined the evolution equations for truncated moments 
of parton distribution, and given their explicit solution to NLO in the nonsinglet 
sector. From a theoretical viewpoint, this fills an obvious gap in the available abun- 
dant literature on QCD evolution equations and the methods for solving them. 
From a phenomenological viewpoint, our results are a useful addition to the set of 
tools available to extract information from the data on scaling violations, and in 
particular provide a new way of dealing with the well-known problem of working 
around our ignorance of the small-x behavior of structure functions. A preliminary 
extraction of a s (M z ) from the scaling violations of the nonsinglet second moment 
looks very promising. While we postpone to future work a fuller phenomenolog- 
ical analysis, we encourage experimental collaborations, which necessarily have a 
much better control of the experimental systematics, to use the simple technique 
presented in this paper as a means to present data on the moments of the gluon 
distribution, as well as to obtain determinations of a s , which would be significantly 
less dependent on model assumptions, in comparison to those obtained with more 
standard techniques. 



Appendix 

We list here a few useful properties of triangular matrices. Consider a generic 
n x n upper triangular matrix T n , with matrix elements a y -, (a y - = for i > j). It 
is straightforward to show that: 

a) The matrices T n form a proper subgroup of GL{n). 
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b) The eigenvalues of T n coincide with the diagonal entries djj. This can be seen 
noting that the secular equation receives contributions only from the diagonal, 
since all other minors of T n — Al vanish. 

c) Let E n be the matrix of right eigenvectors of T n , arranged in columns. E n 
is also an upper triangular matrix, which can be chosen to have all diagonal 
elements equal to unity. Specifically, let (3^ be the j-th right eigenvector of 
T n , defined by 

T n ^) = ajj pU) f (24) 

and let (3^ be the 2-th component of the j-th eigenvector. Then one can 
choose (3^ = 1, and one sees that (3^ = for i > j. Furthermore, the matrix 
E n , with elements E n>i j = f3^\ satisfies 

E~ l T n E n = &i&g(a 3j ) . (25) 



d) The nonvanishing elements of E n satisfy the recursion relation 

= ~ J _— E M?' > ( 26 ) 

a jj a n p=i+l 

initialized by /3j = 1. This recursion relation can actually be solved explicitly 
in terms of minors of the matrix T n , however for our purposes the recursion 
relation itself is more useful, since it can easily be implemented in an evolution 
program. 

e) The inverse matrix E~ l is also upper triangular with unit diagonal entries, 
and can be constructed noting that, for a generic nonsingular square matrix, 
the matrix of right eigenvectors (arranged in columns) is invertible, and the 
inverse is the matrix of left eigenvectors (arranged in rows). This leads to 
a recursion relation for the elements of E' 1 , analogous to Eq. (PBj) . Setting 
(E-% = 7j (j) we find 

7f ) = - i -E7? ) a ra , (27) 

Un Ujj p= i 



It follows that triangular matrices can be diagonalized without computing any 
determinants. The implementation of the recursion relations Eqs. ( p6| - p7| ) is in fact 
so efficient that the entire solution of the evolution equations (for reasonably sized 
matrices, say n < 7 — 8) can be performed analytically, for arbitrary xq, and the 
chosen value of xq substituted only at the end. 
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